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ABSTRACT 


This  report  summarizes  the  results  of  various  analytical  studies  made 
under  Contract  AF  04(611)8018  . 

RELATIONSHIP  OF  EXCESS  BURNING  AREA  AND  CHAMBER  PRESSURE. 
The  relationship  of  chamber  pressure  to  an  arbitrary  increase  in  burning 
area  has  been  derived  for  propellants  of  various  burning  rate  exponents  by 
employing  the  basic  relations  of  solid  rocket  motor  combustion. 

STRESS  ANALYSIS  OF  THE  PROPELLANT  GRAIN.  A  method  of  analysis 
has  been  developed  to  determine  the  elastic  stresses  in  a  propellant  grain 
subject  to  a  constant  body  force  and  uniform  temperature  change.  By  two 
dimensional  analog  of  a  grain,  a  "stress  concentration"  factor  for  the  end 
effects  has  been  determined  to  be  approximately  2-1/4  times  the  plane 
strain  solution.  Numerical  results  are  presented  for  the  120-inch  diameter 
motor. 

STRESS  CONCENTRATION  DUE  TO  A  VOID  IN  THE  PROPELLANT  GRAIN. 
The  effect  a  void  in  the  propellant  grain  has  upon  grain  stresses  has  been 
developed. 

PRESSURE  DISTRIBUTION  IN  A  PROPELLANT  GRAIN  CRACK.  An 
investigation  has  been  made  of  the  pressure  distribution  due  to  burning 
within  a  crack  in  the  propellant  grain.  For  crack  depth/width  ratios 
exceeding  300,  sonic  flow  at  the  crack  exit  is  predicted. 

MOTOR  CASE  HEATING.  The  problem  of  motor  case  heating  due  to 
premature  web  burnout  has  been  investigated.  Predicted  values  of  motor 
case  temperature  versus  time  and  insulation  thickness  are  shown. 

DETERMINATION  OF  SOLID  PROPELLANT  BURNING  AREA  FROM 
PRESSURE— TIME  TRACE.  A  technique  has  been  developed  to  evaluate 
the  results  of  motor  firings  and  calculate  the  actual  burning  area  as  a 
function  of  time. 
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1.0  INTRODUCTION 

In  the  course  of  satisfying  the  technical  requirements  of  Contract 
AF  04(611)8018,  a  number  of  analytical  studies  were  conducted.  The 
purpose  of  this  Technical  Note  is  to  present  the  results  of  the  studies  made 
to  date.  For  the  most  part,  these  studies  have  been  concerned  with  the 
defect-failure  relationships  in  large  solid  propellant  motors. 
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2.  0  TECHNICAL  DISCUSSION 

2.  1  RELATIONSHIP  OF  CHAMBER  PRESSURE  TO  EXCESS  BURNING 
AREA 


The  flow  rate  through  a  choked  rocket  motor  nozzle  is 


w 


CP  A 
c  e 


(1) 


The  flow  rate  from  the  burning  propellant  surface  area  is 


w  =  p,  A  r 
b  s 


and  assuming  the  experimental  burning  rate  law  of  r  =  b  P^  , 
then 


w  =  p,  A  b  P 
b  s  c 


(2) 


For  steady  state  conditions,  the  flow  rates  in  Equations  1  and  2 
must  be  equal.  Thus, 
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Letting  the  subscript  "o"  denote  initial  conditions,  then 


A  graphical  solution  of  Equation  3  for  various  burning  rate  exponents 
is  given  in  Figure  1. 

2.  2  STRESS  ANALYSIS  OF  THE  PROPELLANT  GRAIN 
2.  2.  1  Introduction 

An  area  of  considerable  importance  in  the  design  of  solid  propellant 
grains  is  that  of  three-dimensional  linear  elasticity.  Many  of  the 
complex  problems  in  this  area  are  extremely  involved  and  time 
consuming  or  they  do  not  have  known  exact  solutions.  Therefore, 
simplifying  assumptions  or  simpler  approximate  solutions  are  often 
necessary. 

The  effects  of  the  ends  on  the  stress  distribution  in  a  solid  propellant 
grain  due  to  thermal  and  gravitational  loads  is  one  of  these  complex 
problems.  The  application  of  classical  two-dimensional  elasticity 
in  the  analysis  of  a  grain  acted  upon  by  these  loads  is  well  reported 
(References  1,  2,  and  3).  In  order  that  these  simplified  solutions 
apply,  an  assumption  is  made  that  the  area  to  be  investigated  is  far 
from  the  ends  of  grain.  Thus  a  condition  of  plane  strain  is  assumed. 
But  for  the  most  part,  these  simplified  analyses  do  not  apply  near  the 
ends  of  the  grain.  Since  it  is  well  known  to  propellant  manufacturers 
that  bond  failures  near  the  end  of  a  vertical  standing  propellant  grain 
occur  occasionally,  a  more  detailed  study  is  required  for  this  very 
practical  problem. 
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Some  attempts  have  been  made  to  obtain  three  dimensional  elastic 
solutions  to  the  end  effects  problem.  Knauss  (Reference  4)  suggests 
that  a  variational  method  can  be  used  if  the  form  of  the  radial  and 
axial  deflections  are  assumed.  By  using  the  Fourier  method  suggested 
by  Pickett  (Reference  5),  a  form  for  the  deflections  is  assumed.  The 
behavior  of  the  stress  components  along  the  boundary  are  then  determined. 
Conte,  Miller,  Sensenig  (Reference  6)  and  Messiner  (Reference  7) 
attempted  to  solve  the  end  effects  problem  through  numerical  pro¬ 
cedures.  Difference  equations  replace  the  differential  equations  of 
elasticity  at  discrete  points  in  the  grain.  The  stress  components  are 
predicted  by  using  iterative  techniques. 

None  of  these  suggested  solutions  provide  complete  correspondence 
to  the  original  problem.  If  the  form  of  the  deflections  are  chosen  as 
suggested  by  Knauss,  complete  satisfaction  of  equilibrium  and/or 
the  original  boundary  condition  is  not  attained  (i.  e.  deflection  or 
stresses  along  certain  boundaries  are  not  maintained  at  specified 
values).  Because  of  accuracy  difficulties,  the  iterative  techniques 
also  do  not  completely  satisfy  the  original  boundary  conditions. 

Therefore,  a  better  understanding  of  end  effects  is  required. 

This  presentation  is  separated  into  two  sections.  The  first  part 
deals  primarily  with  the  practical  answers  of  bond  stresses  in  a 
propellant  grain.  The  second  part  presents  the  formal  solution  to 
the  end  effects  problem. 

2.  2.  2  Stresses  in  a  Propellant  Grain 

Propellant,  standing  vertically  being  cured,  can  be  ideally  represented 
as  an  infinitely  long  hollow  cylinder  bonded  at  its  outer  radius  to  an 
elastically  rigid  case.  The  propellant  is  assumed  to  be  elastic. 
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It  is  assumed  that  superposition  principles  remain  valid.  Therefore, 
the  thermal  loading  resulting  from  cure  and  the  gravity  loading  are 
assumed  to  be  separable. 


The  stress  free  condition  for  the  thermal  loading  is  assumed  to  be 
at  the  cure  temperature,  T.  As  a  result  of  the  differences  in  co¬ 
efficients  of  linear  expansion  between  the  propellant  and  its  casing, 
thermal  stresses  result  when  a  uniform  change  in  temperature, 
AT,  is  experienced. 


The  thermal  elastic  stress  equations  of  plane  strain  for  the  propellant 
including  the  effect  of  polymerization  are  found  to  be  (Sections  2.  2.  3 
and  2.  2.  4) : 


where 


A  = 


a  )  AT 
c 


(4) 


Near  the  ends  of  the  cylinder,  these  stresses  are  altered  drastically 
from  the  plane  strain  solution.  To  determine  just  how  radical  a 
change  in  stress  distribution,  the  analogous  two  dimensional  elasticity 
problem  (Section  2.  2.  3)  is  examined. 


II 
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At  the  ends  of  the  cylinder  an  effective  increase  in  maximum  stress 
or  stress  concentration  factor  is  assumed  to  be  the  same  as  for  the 
two  dimensional  solution.  The  maximum  stresses  at  the  bond  surface 
have  a  magnitude 


O’  = 

r 


(5) 


Note  that  the  end  effects  produce  a  shear  stress  not  normally  observed 
by  the  plane  strain  solution. 


The  length  of  the  cylinder  can  also  be  taken  into  account  (Section  2.  2.  3), 
thus  giving  rise  to  higher  stress  concentrations.  For  a  length  to  web 
thickness  ratio  of  3,  the  maximum  stresses  are: 
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A  similar  procedure  for  the  gravitational  loading  results  in  the 
expressions : 


<T 

r 


6 


cr 

z 


T 

rz 


0.  870  p  g  b 


0.430  Pgb 


0,  870  p  g  b 


0.77  p  g  b 


(7) 


The  combined  stresses  due  to  gravity  and  thermal  loading  at  the 
outer  radius  is  simply  the  addition  of  Equations  6  and  7: 


0-  =  +  1.  2  A 
r  — 


(r„  =  1.71 


(T  =2.21 

z 


T  =  1.41 
rz 


c-') 


1  j  ±  0*  870  p  g  b 


('■?) 


A  —  ±  0.  430  p  g  b 
a 


(■>) 


±  0.870  pgb  (1 


±  0.  77  p  gb  11 


(-fl 

{■$) 


(8) 


To  obtain  the  principal  stresses,  the  three  normal  stresses  and  the 
one  shear  stress  are  used.  Since  the  shear  stress,  ,  is  zero, 

one  of  the  components  in  the  tB  plane  can  be  assumed  to  be  a 
principal  stress. 
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Assume  ^"0  ^  ®’i 

O  £1  L 


Figure  2.  Mohr's  Diagram 


From  Mohr's  circle  (Figure  2),  the  principal  normal  stresses, 
ir  and  tr  ,  can  be  determined 


The  maximum  shearing  stress  is  therefore  defined  as 


max  2 

The  use  of  the  maximum  shear  criterion  by  no  means  implies  the 
propellant  material  follows  this  failure  phenomenon;  it  is  to  be  used 
only  as  a  guide.  Further  investigation  as  to  the  mode  of  failure  is 
required  before  a  failure  criterion  is  chosen. 
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2.  2.  2.  1  Application  to  120-Inch  Motor  Grains 
^  The  assumed  mechanical  properties  of  PBAN  propellant  and 
physical  dimensions  for  the  120-inch  grain  are  as  follows: 

E  =  300  psi 

a  =  3.  5  X  10  ^  in./in./“F 

AT  =  100“F 

pg  =  0.063  Ib/cu  in.^ 

b  =  60  in. 

a  =  25  in. 

L  =  60  in. 

AV 

~  =  0.005  .  . 

Substituting  the  above  values  into  Equation  8,  the  stresses  at 
the  bond  are  obtained 


=  10.  44  psi 
=  + 11.  14  psi 
=  +  17.  64  psi 
T  =  1 1.  92  psi 

Therefore,  the  principal  stresses  have  values 


=  28.  34  psi 


^2  ~  "1*0  psi 

and 

T 

max  2 


14.  04  psi 


Note  these  maximum  stresses  occur  at  the  top  of  a  standing  grain 
The  stresses  at  the  bottom  tend  to  cancel  one  another. 
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2.  2.  2.  2  Effect  of  Boot  at  End  of  Grain 

Effective  use  has  been  made  of  a  relief  at  the  end  of  a  propellant 
grain  in  a  form  of  a  boot  or  cut  formed  wholly  in  the  propellant 
liner.  This  boot  allows  a  portion  of  the  grain-case  interface  to 
be  free  during  cure.  Upon  completion  of  cure,  the  void  so  created 
is  filled  with  a  compound  such  as  an  epoxy  resin.  This  procedure 
is  designed  to  keep  the  ends  of  the  grain  free  of  dangerous  cracks. 

If  a  body  is  free  to  move  under  a  uniform  temperature  change, 

no  stresses  are  induced  in  the  body.  Stresses,  proportional  to 

a  ,  are  induced  when  the  body  is  completely  restrained.  The 
est 

boot  creates  a  stress  condition  somewhere  between  these  two 
extremes.  Since  the  modulus  of  most  propellants  is  very  low,  the 
state  of  stress  is  more  likely  to  be  as  in  a  free  body  rather  than 
restrained.  This  also  means  that  a  case  bonded  cylinder  of  a 
length, 

L'  =  L  -  2B 

where  L  =  length  of  original  cylinder 

B  =  length  of  boot  or  cut  (one  end)  , 
can  be  assumed  to  stress  free  edges.  From  the  effect  of  length 
analysis  in  Section  2.2.3,  it  is  noted  that  (within  limits)  as  the 
length  of  the  cylinder  decreases,  the  stress  at  the  end  increases. 
Thus,  although  the  stresses  at  the  extreme  end  of  the  cylinder 
have  decreased,  the  stress  at  the  base  of  the  boot  has  increased. 
Figure  3  shows  other  stresses  at  the  interface  for  L/W  ratios 
of  2  and  3  for  thermal  loading;  Figure  4,  those  for  gravity 
loading. 

There  is  one  relieving  factor  when  the  boot  is  installed.  After 
the  grain  has  reached  ambient  conditions  and  after  cure  has  taken 
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Figure  4.  Stresses  Due  to  Gravitational  Loads  Along  a  Fixed  Edge 

Finite  Length  Cylinder 
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place,  polymerization  has  occurred.  The  filling  of  the  boot  void 
reduces  the  reference  or  datum  temperature  for  that  portion  of 
the  cylinder.  Thus  cylic  temperature  variations  will  be  referenced 
to  the  "ambient  fill  temperature"  and  not  the  cure  temperatures. 
Therefore,  the  stresses  are  greatly  reduced  at  the  ends.  It  should 
be  noted  that  the  maximum  stress  condition  has  moved  to  the  base 
of  the  boot.  A  case  debond  may  occur  as  a  result  of  a  temperature 
variation  at  the  base  of  the  boot,  and  yet  the  ends  remain  crack 
free. 

2.2.3  Stresses  in  a  Hollow  Cylinder 

A  cured  circular  cylindrical  propellant  grain  can  be  idealized  initially 
to  an  elastic,  infinitely  long,  hollow  cylinder  bonded  at  its  outer  surface 
(r  =  b)  to  a  steel  case  (Figure  5).  The  inner  surface  (r  =  a)  is  assumed 


RIGID  CASE 


PROPELLANT 

Figure  5.  Grain  Configuration 
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to  be  free  of  normal  and  shearing  stresses.  At  a  temperature,  T  , 

the  cylinder  is  free  of  stress.  The  cylinder  is  assumed  to  be  elastic 

(having  a  modulus  of  elasticity,  E)  and  obeys  the  Duhamel -Neumann 

(Reference  8)  stress -strain  relations.  Its  coefficient  of  linear 

expansion  is  a  ,  while  that  of  the  steel  case  is  a  .  The  change 
P  ^ 

in  volume  due  to  polymerization  of  the  propellant  can  be  expressed 
A  V 

by  the  ratio  .  The  problem  is  to  determine  the  maximum 
stresses  in  the  cylinder  caused  by  the  decrease  in  temperature  from 
T  to  T  -  A  T  . 


Due  to  a  difference  between  the  coefficients  of  linear  expansion  of 
case  and  propellant,  a  stress  distribution  will  arise  in  the  propellant 
grain  when  the  temperature  changes.  Assuming  the  case  to  be  elas¬ 
tically  rigid,  the  equations  for  the  stress  field  of  an  infinitely  long 
cylinder  can  be  easily  obtained  (Reference  9).  The  three  normal 
stresses  have  the  form  of  the  Lam/ equations.  For  Poisson's  ratio 
1/  =  —  (Section  2.  2.  4)  the  stresses  are 


cr 

r 


(a 


a  )  AT 
c 


(11) 


a 


e 


(q  -  a  )  at  +  — 
pc  3  V 


(12) 


cr 

z 


AV 

(a  -  a  )  AT  +  — 
pc  3  V 


(13) 


If  the  cylinder  is  now  assumed  to  be  of  finite  length,  the  ends  of  the 
cylinder  have  normal  stresses  in  accordance  with  Equation  13.  To 
duplicate  the  stress  condition  of  the  propellant  grain,  no  external 
forces  or  stresses  are  applied  to  the  ends  of  the  hollow  cylinder. 
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Therefore,  the  end  surfaces  must  have  the  normal  stresses  described 
in  Equation  13  removed.  This  is  done  by  superposing  equal  and  opposite 
normal  stresses  on  the  ends  of  the  cylinder.  This  latter  problem, 
sometimes  referred  to  as  the  "reduced  problem,  "  has  no  known  exact 
solution  at  present. 

Therefore,  to  develop  a  more  clear  understanding  to  the  stress  distri¬ 
bution  in  the  cylinder  for  the  reduced  problem,  consider  the  case  of  a 
semi-infinite  slab  (Figure  6).  The  correspondence  between  the  semi¬ 
infinite  slab  and  the  semi-infinitely  long  hollow  cylinder  is  geometric. 
That  is,  if  the  radius  to  the  centroid  of  the  longitudinal  cross  sectional 

area  of  a  hollow  cylinder  becomes  infinitely  large  or  if  the  ratio  of  the 

b 

outer  to  inner  diameters,  —  ,  approaches  one,  the  hollow  cylinder 

degenerates  into  the  slab.  At  the  ratio  —  =  1  ,  the  elastic  stress 

a 

solutions  of  the  two  bodies  would  also  be  the  same. 


R-20<21 


Figure  6.  Semi-Infinite  Slab 
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2.  2.  3.  1  Stresses  in  a  Semi-Infinite  Slab  (Thermal  Loading) 

The  slab  is  assumed  to  have  the  same  elastic  characteristics  as 
the  hollow  cylinder.  If  the  slab  were  imagined  to  be  free,  as  a 
result  of  a  uniform  temperature  change,  AT ,  the  body  would 
have  strains  developed 

€  =€  =€  =qAT 

X  y  z  p 

(14) 

€  =0 
xy 


and  the  steel  case  would  have  strains 


€  =e  =€  =aAT 
X  y  z  c 

(15) 

€  =  0 
xy 

Consider  plane  strain  to  exist  in  the  z  direction.  Since  the  relative 
displacement  between  the  slab  and  steel  casing  must  be  zero,  a 
stress  must  be  developed  throughout  the  body  having  one  magnitude 

0-  =  E  (a  -  a  )  AT  =  Ea  AT  (16) 

X  p  p  c 

As  for  the  cylinder,  the  end  surfaces  must  have  the  normal  stresses 
described  in  Equation  16.  To  remove  these  stresses,  equal  and 
opposite  normal  stresses  are  superposed.  Again  the  "reduced 
problem"  is  defined. 

At  the  end  plane  x  =  0  (Figure  6),  a  uniform  normal  stress  is 
applied 

0-  =  -  E  aAT  .  (17) 

X 
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Note,  a  convenient  form  in  notation  is  introduced 


f  =  £  (x),  f  >  =  ^  ,  £  "  =  , 

n  n  n  ox  n  «  2 


The  second  equilibrium  condition  states 


da-  dT 


=  0  . 


Thus,  the  stress  in  the  y  direction,  ,  can  be  expressed  as 

0  O'  f  '  y 

V  \  n  ^ 


3y  ^  n 
n=l 


m  -  n+1 

y 


(n+l)  ^2^^^^ 


To  satisfy  the  boundary  condition, 


0-  {x,  0)  =  0 

y 


we  obtain 


C^(x)  =  0  . 


To  obtain  a  solution  for  f  ,  the  variational  approach  (Reference  11) 

n 

to  minimize  the  strain  energy  expression  is  used.  The  complementary 
energy  expression  for  plane  strain  has  the  form 


#/7  I 

X,  y,  V 


1 


Z  ^  Z  Zv 

(j-  +  <r  -  - —  0*  O’  + 

X  y  1-1/  X  y 


a  d  f 

— — ‘—I  dxdy  -  I  cr  u« 
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The  second  integral  is  that  portion  of  the  surface  of  the  body,  F, 
over  which  the  displacements  are  specified.  Along  the  clamped 
edge,  the  displacement  is  specified  as  zero.  Thus,  the  integrand 
of  the  second  integral  is  zero.  In  addition,  the  remainder  of  the 
above  expression  can  be  normalized  such  that 


X  = 


X 

h 


(29) 


and 


y 


I 

h 


Introducing  Equations  18,  21,  and  26  into  Equation  28  and  using  the 
variation  principle  for  each  function,  f^  ,  a  set  of  differential 
equations  will  be  obtained  in  the  form 


This  condition.  Equation  30,  is  a  necessary  condition,  insuring 
that  the  integral 


U 


■■/I 


F{f,,  f 


r  ^2’ 


f  )  dy  dx 
n 


minimum 


(31) 


However,  imposing  minimum  complementary  energy  is  not  sufficient 
to  insure  that  the  clamped  edge  (y  =  1)  conditions  are  maintained 
as  necessary  conditions. 
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For  the  displacement  boundary  condition  at  the  fixed  edge,  the 
stress -strain  relations  in  the  x  direction  can  be  expressed  as 

E€  =  (1-V^)  -  0-  1/(1  +  I/)  =  0  .  (32) 

X  ox  x  y 

thus,  a  relation  between  the  normal  inplane  stress  components 
results 


(T 

X 


cr 

y 


(33) 


In  addition,  the  displacement  in  the  y  direction,  v,  is  zero  at 
the  fixed  edge.  Again  using  stress -strain  relationship. 


if 

ay 


a 

X 


1/(1  + 1/) 


(34) 


Substituting  Equations  18  and  26  into  34  and  integrating,  the 
displacement  is  defined  as 

.A  A  f  7""' ■'(!  +  ■') 

''  =  A  ''n(n.l)(n.2)  '  -  *  ''W  ■ 


n=  1 


n=l 


n+  1 


where 


V(x) 


f"  (I  -  v^) 
_n _ ' 

n  (n+l)(n+2) 


f  V(  \  +  V) 


n+  1 


Note  that  the  displacement  boundary  condition  in  the  y  direction 
is  satisfied  identically  by  the  function  V  (x)  .  The  displacement 
in  the  x  direction,  u  ,  or  its  derivative,  however,  as  described 
by  Equation  32,  is  not  assured.  Thus,  the  relationship  in  Equation  33 
must  be  used  as  an  auxiliary  boundary  condition  (Reference  12). 
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Substituting  Equations  18  and  26  into  Equation  33,  the  non-holonomic 
condition  is  provided  at  the  fixed  boundary  (y  =  1) 

(36) 


111  111 

E  £  -  E 

n  l-v  n{n+l) 


n=l  “  ‘  '  n=l 

or  simply  expressed  in  functional  form 


G(f  ,f'')  =  0  . 

n  n 


(37) 


Now  to  insure  that  Equation  33  is  satisfied  as  well  as  the  comple¬ 
mentary  energy  being  a  minimum,  the  Lagrangian  multiplier,  ^  , 
(Reference  12)  is  used  in  the  variation.  Thus,  Equation  31  is 
modified,  such  that, 

U*  =  U  +  |G  .  (38) 

The  variation  of  the  modified  energy  can  be  expressed  as 


where 


6  U*  =  6  (U  +  |G) 


^  =  l(x) 


(39) 


U 


The  complementary  energy  can  be  determined  by  substituting 
Equations  18,  21,  and  26  into  28 


/ 


m  m  „  „_n+i+2 

^  ^  -n+i+2  ^  y  y  Vn^ 

i  n^  in(n+l)(i+l) 

1-1  n=  1 


m 


m  £  i£ «  ^+i 
i=l  n=l 


I  n 


21/ 

\-v 


f.f'v 
1  n 


n+i 


i=l  n=l 


n(n+l) 


(40) 


dy  dJT 
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However,  to  introduce  the  subsidiary  boundary  conditions  according 
to  Equation  38,  the  following  relation  must  hold 


U*  =  U  + 


II  'it  ^n-7^  §^)| 


Since  the  functions  f  ,  f.,  and  i  are  functions  in  x  only,  the  y 

n  1  ’  ^ 

variable  can  be  evaluated  for  y  ^  =  0  and  y^  =  1 


m  m 


II  c  II 


L  I  iMn.lKHDIn+W) 


m  m 


2 


•  (  • 


f.'  f 


m  m 


+  -=-  /  Y  y  V 

^  ^  in{n+i+l)  \-v  ^ 


f.  f  " 

1  n 


l-V  4 


i=l  n=l 


i=l  n=  1 


n(n+l){n+i+l) 


m 


m 


I  V  f  .  liL  V  Ja 

^  4^  n  \-v  4-i  n(n-t 
1=1  1=1 


f  " 
n 

n{n+l) 


dx 


(42) 


Applying  Equation  30  to  42  for  each  function  f.  and  |  we  obtain: 


6  U*  = 


4 

1  -V 


m  m 

1=1  n=l 


+  V 


i=l 


i(i+l)  +  (n*l)  n  \  „ 

in(n+l)(i+l)(i+n+l)  /  n 


m  m 


i=l  n=l 


f.  dx  =  0 
1 


_ 2 _ 

i  n(n+l)(i+l)(n+i+3) 


(43) 


n 
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It  is  noted  that  the  function  within  the  braces  has  the  general  form 
of  the  Euler-Lagrangian  equation. 

a.  +  b.  f."  +  c.  f.^  =  b  (44) 

in  in  1  in  1  | 

There  are,  however,  m  number  of  Euler  equations  that  must  be 
solved  simultaneously. 


fa..  +  b..  +  c..D^1f.  =  0  (45) 

L  ij  J  1 

where 

i,j  =  (1.2, 


Solution  to  the  previous  equation  is  obtained  by  merely  replacing  the 
derivative  with  its  characteristic  root  or  eigenvalue,  X.  _  ,  and  solving 
the  resultant  determinant  (D  matrix).  (Refer  to  Reference  13. )  The 
form  of  the  functions  is  therefore 

.X 

e  ^  .  ( 46 ) 

As  a  numerical  example,  assume  i  =  3;  the  resultant  D  matrix 
would  be  expressed  as 


•  •  •  ,  m) 

,  etc. 


2 

(i  +  .j)  ^ 

I 

(l  +  if' 

\  ^ 

2 

(± ,  + 

2^  , 

'1 

-  [j  *  '7  ~ 

10  ' 

1 

(26 

)  i-r 

36 

3 

ys  2/1-1' 

84  ’ 

2(1-1') 

1 

/I  +  i!r')  ^ 

2 

0“ 

1 

(i  +  iif)  ^ 

'2 

(z  6/1-1' 

36  ' 

3 

v  3 ; 

1  -r 

126 

’  2 

\9  72  /  1  -  v 

288 

■  ’  6(1-1/)  ° 

2 

.  (±.  +  Ajl. 

+  2.^ 

1 

.  (i  +  ii!? 

1  dL  t 

2 

[±  t  + 

^3 

T 

Vl5  i/i-v 

84  ' 

2 

^9  *  72  , 

1  1 -If 

288 

•  5  ■ 

\63  9/  1  -v 

648  ' 

12(1-1/)  ° 

1 

_ t 

n*  1 

V 

-.2 

i 

2(1-1')  ' 

mT 

D  «  1  * 

12(1-1') 

L/  t 

0 

=  0 


(47) 
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Thus,  the  form  of  the  functions,  f.  and  i,  can  be  expressed  as 


6 


X  .X 
1 


f.  = 


e  = 


X  .X 


B  .  e 
li 


(48) 


E 

i=i 

6 

E 

i=l 


""li' 


X.  X 


For  convenience,  a  Poisson's  ratio  of  l/2  will  be  assumed  as  a 
sample  of  this  procedure  since  it  is  representative  for  most  solid 
propellants.  The  eigenvalues  are  listed  in  Table  I. 


TABLE  I 

EIGENVALUES  OF  EQUATION  47 


=  -0.69357712 

=  -1.7043982  + 

i  2.849476 

^4,5 

=  -3.5030151  + 

i  6. 7263859 

=  -14.030111 

Since  the  eigenvalue,  X  ,  is  a  root  of  the  determinant,  so  is  its 
conjugate,  X,  as  well  as  their  negatives.  Since  Equation  46  is 
of  the  exponential  form,  a  decaying  or  negative  eigenvalue  will 
therefore  be  used.  Thus,  the  upper  limit  of  the  summation  on 
f.  is  twice  the  number  of  functions  chosen. 
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By  substituting  Equations  48  into  Equation  47,  moving  the  first 
column  matrix  of  Equation  45,  A. . ,  to  the  right  hand  side  of  the 
equation,  and  using  n-1  ,  homogeneous  equations,  the  ratio  for 
each  element  in  the  matrix  (B/A).,  and  (d/A)  is  obtained 


’b/a' 

■=11 

=11 

°ir 

-1 

c/a 

= 

®21 

^21 

°21 

• 

^21 

d/a 

1 

,®31 

s. 

°31_ 

/31 

The  axial  stress  at  x  -  0  is  normalized  so  that  v  =  a  ET  =  1 

and  T  =  0  .  * 

xy 


In  satisfying  the  remaining  boundary  conditions,  therefore,  ^he 
following  equations  must  be  satisfied: 
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or  substituting  values  from  Equation  49 


1 

1 

1 

1 

1 

1 

■^r 

’f 

(B/A)j 

(B/A)^ 

(B/A)j 

(B/A)^ 

(B/A)^ 

(B/A)^ 

^2 

0 

(C/A)| 

(C/A)^ 

(C/A)j 

(C/A)^ 

{C/A)^ 

(C/A)^ 

0 

S 

^3 

^4 

^5 

= 

0 

A  3 

^  4 

^  5 

■^6 

0 

1 

^  3 

A  4 

^  5 

>JP 

O' 

1 

0 

The  values  of  the  coefficients  are  summarized  in  Table  II. 

TABLE  II 

SOLUTION  VECTOR  OF  EQUATION  51 
Aj  =  1.4737285 

A  =  -0.60773211  +  i  0.45652196 

A.  _  =  0.60715210  +  i  0.44805945 

4, 5  - 

A^  =-0.47256973 

By  substituting  the  coefficients  from  Table  II  into  Equation  48  and 
into  Equations  18,  21,  and  26,  the  stresses  along  the  fixed  edge 
are  obtained  (Figure  7).  To  obtain  the  final  stress  distribution 
along  the  wall.  Equation  16  must  be  added  to  Equation  18.  This 
is  represented  in  Figure  8. 
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Figure  7.  Stresses  Along  a  Fixed  Edge  for  Reduced  Proble 
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Figure  8.  Stresses  Along  a  Fixed  Edge  for  Thermal  Loading 
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2.  2.  3.  2  Stresses  in  a  Semi -Infinite  Slab  (Body  Forces) 

The  second  consideration  is  that  of  the  stresses  induced  by 
gravitational  body  forces.  The  stresses  in  a  vertical -standing, 
hollow,  infinitely  long  cylinder,  as  in  Figure  5,  can  again  be 
easily  determined  (Reference  1). 

T  =  ^  r  fl  .  (if 
rz  2  L  O 

a  =0-  =0-^  =  0 

r  z  0 

All  boundary  conditions  are  satisfied  except  for  the  ends.  Again, 
the  plane  strain  solution  of  a  semi -infinite  slab  will  be  considered. 
The  expression  for  the  maximum  shearing  stresses  can  be  readily 
found  to  have  the  form 

=  pgy  •  (53) 

From  Equation  21,  it  is  obvious  that  the  same  eigenvalues  of 
Equation  47  can  be  used  and  only  Equations  50, and  51  will  change. 

To  satisfy  the  boundary  condition,  the  following  equations  must  be 

satisfied;  this  time  the  shear  stress  T  =  pgh  =  1  is  normalized 

xy  ® 

at  the  fixed  boundary  (y  =  1) 


6  6 


6  6 


6  6 
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Substituting  values,  as  related  in  Equation  52,  into  Equation  54: 


1 

1 

1 

1 

1 

1 

'V 

”  0~ 

(B/A)^ 

(B/A)^ 

(B/A)^ 

(B/A). 

(B/A)^ 

^2 

0 

(C/A)j 

(C/A)^ 

(C/Alj 

(C/A)^ 

(C/Ajj 

(C/A)^ 

0 

X 

4 

= 

1 

*  2 

A  3 

*  4 

*  5 

(vf) 

■^6 

0 

^  1 

"  2 

^3 

'’^A> 

*4 

*5 

*6 

0 

The  values  of  the  coefficients  are  summarized  in  Table  III. 

TABLE  III 

SOLUTION  VECTOR  OF  EQUATION  55 

Aj  =  2.7651613 

A^  3  =  0.43679095  +  i  0.70784614 

A^  3  =  0.57126581  +  i  0.59331908 

A^  =  -0. 54546726 


Substituting  the  coefficients  from  Table  111  into  Equations  18,  21, 
and  26,  the  stresses  along  the  fixed  edge  can  be  obtained, 

(Figure  9).  Adding  Equation  53  to  the  results  from  Figure  9,  the 
final  stresses  along  the  fixed  edge  can  be  found,  as  represented 
by  Figure  10. 
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Figure  9.  Stresses  Along  a  Fixed  Fdge  for  Reduced  Problem 


2.  2.  3. 3  Stresses  at  the  Re-Entrant  Corner 
As  an  alternate  check,  the  solution  for  the  stresses  at  the  point 
(r,  z)  =  (b,  o),  for  the  hollow  cylinder,  may  be  obtained  directly 
without  solving  the  complete  problem  and  without  resorting  to  the 
variational  method  just  described.  The  pertinent  boundary  conditions 
for  the  cylinder  under  thermal  loading  are 

O’  (a,  o)  =  0 
r 

<r  (r,o)  =  0  (56) 

z 

T  {r,o)  =  0 
rz 

By  the  Duhamel- Neumann  stress-strain  relations 


(T  r  i/(a  +  0-  )  -  E  a  AT  +  E  € 
a  z  r  P  P  P  “ 


(T  =  V{<T  +  or  )  -  E  a  AT  +  E  € 
z  TO  P  P  p  z 


Eliminating  o-  from  the  above  equation,  cr  is  obtained 
r  “ 


(57) 


.  ^ 

■'e  =  *  1 


(58) 


Eliminating  o-^  from  Eqatuion  57,  is  obtained: 


,  „  E  a  AT  E 

0-  =  -  <r  +  — 2 - - 2  (€«!'  +  €  )  .  (59) 

r  V  z  V  i/(  1  + 1/)  6  z' 


Since  the  cylinder  will  experience  equal  strains  due  to  the  thermal 
expansion  o^  the  elastically  rigid  boundary  and  can  be  expressed 
as 

€  =  €  =  ,  (60) 

V  z  c 
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Equation  59  can  be  rewritten  in  the  form 


(T 

r 


1  -V 
V  Z 


E 

V 


a  )  AT 
c 


(61) 


The  stresses  at  the  end  of  the  cylinder  can  then  be  evaluated  by 
introducing  Equations  56  into  61,* 


and 


a 

r 


E  (a  -  Q  )  AT 
P  P  c 

V 


(62)  ; 
(63) 


If  the  plane  strain  solution  is  considered  for  ^  ~  J  •  f^quation  61 

I 

can  be  written  as 


a  =  0-  +  2E  (a  -  Q  )  AT  .  (64) 

r  z  p  p  c 

AV 

From  Equation  13,  neglecting  .  the  radial  is  expressed  as 

(T  =  -  E  (q  -  a  )  at  fl  +  (-)^1  +  2E  (a  -  a  )  AT 
r  ppc  i  a.  J  ppc 

2 

=  -  E  (a  -  a  )  AT  [(-)  -  l] 

ppc  La  J 

which  corresponds  to  Equation  1 1  evaluated  at  r  =  b  . 


From  Equation  34  of  the  semi -infinite  slab  solution  and  Equation  65 

of  the  cylindrical,  reduced  problem,  i.  e. ,  o-  (b,  o),  is  expressed 

z 

by  Equation  13,  the  two  equations  are  of  the  same  form.  The 
radial  component,  r  ,  corresponds  to  the  y  direction,  and  the 
axial  component,  z  ,  corresponds  to  the  x  direction  of  the 
earlier  problem.  Thus,  at  the  one  point  (r,z)  =  (b,  o)  or 
(x,  y)  =  (0,  1),  the  solution  considering  plane  strain  corresponds 
exactly  to  the  cucially  symmetric  solution. 
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The  remaining  question  to  answer  is  whether  the  plane -strain 
solution  for  other  points  in  the  semi-infinite  slab  can  be  an 
analogous  solution  for  the  cylinder.  Goodier  (Reference  14) 

pointed  out  that,  within  the  neighborhood  of  a  boundary,  the  field 

• 

equations  of  elasticity  in  cylindrical  coordinates  for  axially 
symmetrical  stress  systems  reduce  to  the  field  equations  of  plane 
strain.  (A  separate  proof  of  this  is  presented  in  Section  2. 2.  5  . ) 
When  all  values  of  r,  for  the  entire  body,  are  large,  this  state¬ 
ment  can  be  extended  for  the  entire  section.  The  stress  distribution 
is  analogous  to  the  plane  strain  solution,  provided  the  same 
boundary  conditions  exist.  Therefore,  this  analogy  can  be  extended 
to  a  cylinder  of  finite  length,  b,  and  the  corresponding  solution 
of  an  infinite  strip  of  length,  b  . 


2.  2.  3.  4  Effect  of  Length 

The  length  of  the  infinite  strip  can  be  taken  into  account  by  simply 
replacing  L  in  Equation  46  with  assumed  expressions 


f. 

1 


f 


2 


3 


A.'  cosh  X. .  X 
1  1 


B.*  cosh  \  X 

X  1 


C.'  cosh  X. ,  X 
1  1 


(66) 


I 


D.'  cosh  X  .  X 
1  1 


This  form  can  be  assumed  since  it  is  also  a  solution  to  the  Euler- 
Lagrangian  equation  as  discussed  in  Section  2.  2.  3.  1.  Thus,  the 
reference  coordinate  system  is  changed  to  the  center  of  the  body. 
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Substituting  Equations  66  into  18,  the  stresses  in  the  x  direction 
can  be  determined  and  referenced  to  the  center  of  the  body.  The 
same  can  also  be  done  for  the  shear  stresses,  using  Equation  21. 
To  satisfy  the  boundary  conditions,  the  following  must  be  true; 


COih  \  j  - 


i.~)  coih  X  j 
A  ]  12 


(S,^co.hX,| 


X  j  linh  X  j  j 


iinhX^I 


•‘‘“‘Si 


co»h  X ,  - 
1  2 


(^)  co«h  X .  ^ 
A  j  1  2 


(^)  coih  X  .  ^ 
A  i  I  2 


X  .  *inh  X .  — 
1  >  2 


X  (^)  linh  X  .  ^ 
i  A,  i2 


X  (^1  *inh  X .  j 
i  A  .  ,2 


b  — 1 


co»h  X  — 
m  2 


(^)  co«h  X  ^ 

Am  . 


A  co.hX  ^ 
Am  m2 


X.  iinh  X  •- 
m  m2 


X  (f)  linh  X  ^ 
m  A  ^  m2 
m 


^  (T> 

m  A  ^  m2 

m 


J 


1 

^2 

0 

A. 

1 

= 

* 

A 

0 

m 

(67) 


Thus,  substituting  the  above  Equation  67  into  Equations  18,  21,  and 
26,  for  b  =  1  and  v  =  ^  ,  the  stresses  are  obtained  along  the 
fixed  edge  (Figure  11).  Employing  the  same  methods  for  the 
gravity  loading,  the  stresses  along  the  fixed  edges  may  be  obtained 
(Figure  12). 


2.  2.3.5  Discussion 

Although  the  plane  strain  solution  is  basically  two  dimensional  for 
the  slab,  the  results  are  analogous  to  the  hollow  cylinder.  Main¬ 
taining  plane  strain  throughout  the  semi-infinite  slab  assures  that 
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a  relationship  exists  among  the  stress  components.  Since 
compatibility  conditions  must  hold,  a  stress  reversal  of  the  normal 
stress,  O'  ,  at  the  end  is  required.  The  form  of  the  stresses 

y 


A  e  •  f  (y) 

n  n 


The  form  is  similar  to  that  predicted  by  Horvay  (Reference  15) 
and  Kent  (Reference  16).  Two  experimental  investigations,  Durelli 
and  Tsao  (Reference  17),  and  Zandman,  Post  and  Redner  (Reference 
18),  also  show  that  the  experimentally  derived  stresses  can  be 
represented  by  the  same  form  as  Equation  68. 

The  comparison  of  the  predicted  results  to  Zandman’s  experimental 
results  cannot  be  applied  directly,  since  the  Poisson's  ratio  of  the 
birefringent  plastic  was  not  reported.  The  principal  stress  difference 
at  the  corner  is  reported  to  be 


2  k=o 


From  previous  analysis,  Equation  33,  the  relationship  of  the 
principal  stresses  is  satisfied  only  if  =  -^  .  In  order  to  make 
a  comparison,  it  is  assumed  that  the  plastic  material  has  a 
Poisson's  ratio,  v ,  of  t  • 


The  principal  stress  difference  function  (Figure  13)  follows  quite 
closely  to  the  experimental  results  obtained  by  Zandman.  However, 
the  principal  direction,  9  ,  is  off  as  much  as  20  degrees.  Some 
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of  the  principal  factors  that  may  alter  the  magnitudes  of  the 
direction  function  are: 

1.  The  actual  displacement  boundary  condition  at  the  inter¬ 
face  is  not  rigorously  held.  The  principal  cause  being  that 
the  bonding  agent  is  not  infinitely  rigid. 

2.  It  is  difficult  to  measure  the  isoclinics  near  the  corner 
and  along  a  fixed  edge. 

3.  Relaxation  of  the  stress  within  the  material  due  to  visco¬ 
elastic  properties  causes  a  difference  in  stress  magnitudes. 

4.  The  function  f(y)  may  have  a  greater  number  of  terms 
(i.  e.  ,  order  of  the  function  being  greater  than  three). 

5.  Poisson's  ratio  of  the  plastic  is  other  than  ^  . 

From  the  theoretical  standpoint,  Horvay  (Reference  15)  and 
Williams  (Reference  19)  present  the  argument  that  a  stress  singu¬ 
larity  exists  at  the  re-entrant  corner  and  therefore  the  shear 
stresses  are  infinite.  Mesner  (Reference  7)  agrees  with  this 
argument.  For  two  of  the  arguments  (References  15  and  19),  an 
initial  assumption  of  the  form  of  stress  function  causes  the  stress 
singularity.  The  above  procedures  appear  to  be  correct,  but  do 
not  agree  with  experimental  fact.  From  numerical  studies 
(References  5  and  7),  the  same  prediction  of  high  stresses  results. 
Near  the  re-entrant  corner,  the  numerical  solution  to  the  biharmonic 
equation,  with  mixed  boundary  conditions,  becomes  very  sensitive 
to  small  changes.  Thus,  numerical  accuracies  in  the  method  of 
solution  may  not  provide  adequate  bounds  on  the  iteration  pro¬ 
cedures.  This,  in  turn,  reflects  itself  throughout  the  remainder 
of  the  body,  causing  other  boundary  conditions  not  to  be  satisfied. 
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^.2.3.6  Conclusions 

The  following  conclusions  may  be  drawn  from  the  investigation 
presented  in  Section  2.  2.3: 

1.  The  end  effects  problem  in  a  hollow  cylinder,  at  present, 
does  not  have  a  know  exact  solution. 

2.  The  semi-infinite  slab  analog  of  the  end  effects  problem 
has  a  solution. 

3.  The  variational  energy  methods  alone  do  not  provide  a 
strong  enough  condition  to  insure  that  all  the  boundary 
conditions  are  met.  Boundary  conditions  are  met,  however, 
when  the  Lagrangian  multiplier  is  used  in  conjunction  with 
complementary  energy  methods. 


2.2.4.  The  Lame  Solution 

For  an  infinitely  long,  hollow  cylinder  bonded  to  an  elastically  rigid 
case,  the  stress-strain  relations  are 


I 
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\  ~  V  -  V 
-V  1-1/ 


(T 
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-V  ~v  1 
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AT  + - 

3  V 
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(70) 


and  have  a  value  at  the  bond  radius 


(r=b)  interface  of 


a  AT  + 


1  ^ 
3  V 


From  equilibrium  conditions  with  no  body  forces 


(71) 


8r 


and  knowing  that 


€ 

r 


8u 

9r 


(72) 


(73) 
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and 


(74) 


with  the  added  restriction  that 


the  stresses  for  v  =  —  ,  are  then  found  to  be 


and 


where 


(T  = 

r 


_ r 

2 


(75) 


O',  = 


_ r 

2 

b 


(76) 


P  =  E 
o 


From  Equations  70,  75,  and  76, 


(77) 
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2.2.5  Coordinate  Transformation  Correspondence 

The  equations  of  elastic  equilibrium  of  the  cylinder,  with  no  body 

forces  may  be  written  as 


8e 

9r 
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(1  -  21/) 
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1  9u 
r  9r 
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9e 

9r 
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(1  -  2v) 
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where 
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For  the  coordinate  transformation,  the  first  equation  could  be 
rewritten  as 


8  8u 


8w 

8ti 


u  € 

1  +  T16 


)  +  (1-21') 


4. 

81 


(i  +  UfJ 


0 


Assuming  c— >  8  and,  thus,  €-40  and  neglecting  higher  order  terms 
gives 


8  ■8u 

84  'ae 


8w . 

aV 


+  (1-21/) 


0 


which  is  identical  to  the  plane-strains  solution.  Therefore,  the 
displacements  and  stress  in  a  small  region  of  a  cylindrical  coordinate 
system  are  the  same  as  plane  strain  problems,  provided  the  same 
boundary  conditions  exist. 
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2.3  DEFECT  ANALYSIS 

The  analysis  of  stress  and  strain  in  linear  viscoelastic  bodies  is  considered 
applicable  when  an  analysis  of  a  solid  rocket  motor  propellant  grain  is  in 
question.  A  typical  case  where  a  viscoelastic  solution  can  be  applied  is 
that  of  a  void  defect  in  a  propellant  grain.  As  the  internal  surface  burns, 
the  void  can  only  experience  the  stresses  resulting  from  the  pressure 
loading.  But  suppose,  by  some  means,  the  void  was  to  ignite.  The  occurrence 
of  ignition  prior  to  the  inner  burning  surface  reaching  the  void  would  cause 
added  stresses  in  the  immediately  surrounding  area. 

It  is  assumed,  therefore,  that  the  defect  is  represented  by  a  spherical  cavity 
in  an  infinite  linear  viscoelastic  media.  At  some  time  the  surface  of  the  void 
is  subjected  to  pressure  transient.  Inertia  effects  are  to  be  included  in  the 
analysis.  No  thermal  effects  are  to  be  included  either  in  mechanical  and 
physical  properties  or  in  the  thermal  stress  field.  This  investigation  is  of 
interest  to  support  the  continuing  defect  analysis  under  the  Nondestructive 
Testing  (NDT)  program.  In  addition,  this  presentation  may  show  possible 
areas  of  analysis  not  considered  heretofore. 

2.  3.  1  Stress  Equilibrium 

The  equilibrium  equation  expressed  in  spherical  coordinates  for 

axisymmetrical  loading  condition  has  the  form  (Reference  20): 


if  ^  T  ^  °  ■ 


(78) 


where  R  is  a  body  force. 
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To  include  the  inertia  effects,  the  body  force  can  be  expressed  as 

R  =  -p  .  (79) 

8 

t 

where  p  =  mass  density 

u  =  radial  displacement  . 

In  order  that  a  relationship  between  the  displacements  and  the  stresses 
be  defined,  it  is  required  that  the  special  distribution  of  the  space  and 
time  variations  be  separable.  A  convenient  way  to  remove  this  de¬ 
pendence  is  to  apply  the  Laplace  transform: 

CO 

F  {X,  b)  =  J  f  {X,  t)  e’®^  dt  ,  (80) 

0 

where  f  (X,  t)  is  real  and 

CD 

f(X,  t)  =  ^  J  f  {X,  s)e®^ds  .  (81) 

o 

It  can  be  shown  that  the  removal  of  the  time  variable  allows  the 
application  of  the  already  developed  theories  of  elasticity.  Thus, 
the  extensive  literature  of  elasticity  can  be  utilized  in  the  viscoelastic 
analysis.  Applying  the  Laplace  transform  to  Equation  78  gives 

2  - 

p  8  u  ,  (82) 
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where  (  )  refers  to  the  transformed  state 
and  s  =  operator  . 


From  elasticity  theories  (Reference  5) ,  the  displacement,  u  ,  is 
related  to  the  stresses,  c  ,  as  follows: 


Substituting  Equations  83  into  82,  gives  the  relation 


8^u  28u  2u  p  2  — 

^  X  +  2|ji 


(84) 


Rearranging  terms,  the  well-known  Bessel's  equation  is  obtained 
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where 
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The  solution  to  Equation  85  is  known  (Reference  21)  as 
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Since  the  media  is  assumed  to  be  infinite,  only  the  decaying  function 
will  be  used 


“  c  ,  ..  -asr 

u  =  — y  (asr  +  l)e 

r 


(87) 


The  stresses  can  be  determined  by  substituting  Equation  87  into  83 


r  3 
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—  —  2  2  2 
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To  evaluate  the  constant  c  ,  the  actual  stress  at  the  surface  at 
time  t  =  0  can  be  represented  as 


(^»  \  (Tq.  0)  • 


(89) 


Likewise,  its  transformed  function  can  be  defined  as 


(^»  ®)  =  (^0’  ' 


(90) 


Therefore,  the  coefficient  has  an  expression 
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The  actual  stresses  in  the  body  can  be  expressed  therefore  as: 


_  2  2 
4  p  (1  +  asr)  -  p  s  r 


_  2  2 
4p  (1  +  asr^)  -  ps  r^ 


asrn  +  St 
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2 p  (I  +  asr)  +  X  (asr) 

_  2  2 
4p  (1  +  asr^)  -  ps  r^ 


asr^  +  st 

e  ds  .  (92) 


The  evaluation  of  Equa.tions  92  can  be  performed  by  either  direct  or 
numerical  integration.  However,  in  order  to  evaluate  the  stresses 
the  transformed  Lame  constants  must  be  known. 


For  most  materials,  especially  a  solid  rocket  motor  propellant 
grain,  the  viscoelastic  properties  vary  between  type  or  class  of 
material  so  that  it  is  necessary  to  determine  experimentally  the 
properties  for  each.  In  order  to  apply  any  experimental  data  to 
a  viscoelastic  stress  analysis,  the  response  or  time  interval  of  the 
experiment  must  be  the  same  as  the  interval  of  interest. 

In  addition,  the  type  of  experiment  and,  thus  ,  the  phenomenological 
model  assumed  is  of  interest.  It  has  been  stated  (Reference  22) 
with  good  basis  that  no  solid  propellant  can  be  simulated  by  a  single 
Maxwell  element.  It  would  appear  that  a  number  of  Maxwell  elements 
in  parallel  can  more  adequately  represent  the  dynamic  behavior  of 
propellant  grain.  With  the  additional  as  sumption  of  linear  viscoelastic 
behavior,  an  approximate  model  of  a  finite  number  of  Maxwell 
elements  can  be  used  over  a  given  frequency  spectrum. 
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Welch  (Reference  23)  describes  an  asymptotic  procedure  to  approxi¬ 
mate  the  complex  modules.  This  procedure  is  similar  to  the  Bode 
(Reference  24)  technique  used  in  electrical  network  synthesis. 
Citerley  (Reference  25)  describes  a  similar  procedure,  except  a 
digital  procedure  is  used  in  determining  the  stiffness  of  each  element. 
By  arbitrarily  selecting  the  time  constants  of  the  various  elements  for 
the  particular  range,  adequate  representation  is  developed.  Once  the 
complex  modulus  has  been  determined  for  tensile  as  well  as  shear,  a 
transform  of  the  Lame  constants  is  possible. 

After  the  transformation  is  accomplished  and  the  stresses  in 
Equation  92  have  been  evaluated,  it  remains  to  combine  these  results 
with  a  criterion  for  determining  if  fracture  will  occur.  Unfortunately, 
no  clear  failure  criterion  exists  for  the  propellant  grain  material. 

2.  3.  2  Fracture  Criteria  of  Viscoelastic  Materials 
Fracture  criteria  applicable  to  triaxial  states  of  stress  are  seldom 
available  for  any  material.  Determination  of  such  criterion  for  a 
given  material  is  a  formidable  experimental  task  even  for  elastic 
behavior,  let  alone  viscoelastic.  Uniaxial  data,  however,  is  usually 
abundant.  The  possible  generalizations  of  uniaxial  data  to  multiaxial 
states  are  numerous  and  the  selection  of  a  particular  criterion  is  more 
or  less  arbitrary.  In  this  light,  we  see  that  indiscriminant  application 
of  such  criterion  is  unreliable  and  could  label  as  safe,  stress  states 
that  actually  cause  fracture. 

As  an  alternative  to  trying  to  establish  a  precise  criterion,  upper 
and  lower  bounds  may  be  found  by  using  fundamental  properties  of 
the  fracture  surface  and  available  uniaxial  data. 
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The  fracture  envelope  for  uniaxial  tension  for  various  strain  rates 
at  constant  temperature  appears  as  follows. 


For  a  given  temperature  and  strain  rate,  then,  let  the  tensile 

fracture  strength  be  cr^  ,  Let  the  compressive  strength  be 

cr  .  If  cr  is  not  available,  we  must  assume  cr  =  o-  .  If 
0  0  10 

the  shear  strength  is  given,  for  which  the  principal  stresses  are 
(cr  ,  -<r  ) ,  an  alternate,  but  similar,  procedure  to  what  follows 
must  be  used  to  obtain  the  fracture  spectrum. 

Given  o-  and  v  ,  consider  the  fracture  criteria  to  be  representable 
1  o 

in  a  three-dimensional  space  in  which  the  coordinate  axes  are  principal 
stresses.  The  intersection  of  the  fracture  spectrum!  with  any  coordinate 
plane  appears  as  in  Figure  14,  where  construction  lines  have  been 
included  for  clarity. 

The  assumptions  employed  in  Figure  14  are  that  the  material  is 
isotropic  and  that  the  fracture  locus  varies  linearly  with  hydrostatic 
pressure.  Thus,  the  spectrum  is  only  partial  with  respect  to 
hydrostatic  stress,  but  is  complete  in  the  octahedral  shear  plane. 
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2.  4  DETERMINATION  OF  THE  PRESSURE  DISTRIBUTION  IN  A 
PROPELLANT  GRAIN  CRACK 

2.  4.  1  Analysis 

The  seriousness  of  a  grain  defect  may  depend  on  whether  or  not  the 
defect  will  propagate  or  tear  when  it  is  exposed  to  the  hot  combustion 
gases  and  starts  to  burn.  Before  determining  from  a  stress  analysis, 
whether  the  defect  will  propagate  or  tear  the  pressure  distribution 
along  the  propellant  defect  must  be  determined.  The  purpose  of  the 
following  analysis  is  to  determine  this  pressure  distribution.  The 
analysis  is  applicable  for  a  propellant  grain  crack  with  variable 
cross-sectional  area. 


In  the  analysis,  the  following  simplifying  assumptions  are  used: 

1.  Mass  is  added  at  a  constant  flame  temperature  and 
at  zero  velocity 

2.  The  perfect  gas  equation  of  state  is  considered  valid 

3.  The  divergence  of  convergence  angle  of  the  propellant 
grain  crack  is  considered  small 

4.  The  flow  in  the  propellant  crack  is  quasi  one -dimensional  . 


Considering  a  flow  volume  element  as  shown  in  Figure  15,  the 
continuity  equation  can  be  expressed  as 


(93) 


R-20822 
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the  momentum  equation  as 


^  (wV)  .  ^  (PA)-P^ 


- 


and  the  conservation  of  energy  equation  as 


.  ,  d  ^  m  .  w 

w  h  =  — —  jw  C  T  +  - r — 

p  p  dx  \  p  2 


By  substituting  w  =  p  V  A  into  Equation  93  and  differentiating 
the  product,  then 


.  dV  dA 

dp_  p  '  dx  dx 

dx  ^  VA 


Likewise,  differentiating  the  products  in.  Equation  94  and  solving  for 

^  ^  I  tti6n 


By  using  the  perfect  gas  equation  of  state  and  substituting  Equation  97 

d  P  ,  d  V 

for  A  — —  ,  then  Equation  95  yields,  when  solved  for  — — 

OX  ox 


,,,  (v  -  1)  w  h  +  W  V^  +  yTA  V-yVp4^ 

d V  ^  P  P  \  2  /  p _ S  dx 

dx  (yPA  -  wV) 

(96 

Equations  96,  97,  and  98  were  solved  on  the  IBM-1620.  Since  the 
two  boundary  conditions  are  the  velocity  at  x  =  0  and  the  static 
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pressure  at  x  =  L  ,  an  iterative  procedure  was  used  to  solve  the 
set  of  differential  equations. 

A  typical  solution  of  the  equations  is  shown  in  Figure  16.  In  this 
solution,  the  skin  friction  force  was  neglected  and  the  local 
propellant  burning  rate  was  determined  from  r  =  b  ,  where 
P  is  the  local  static  pressure. 
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Figure  16.  Pressure  Distribution  in  a  Grain  Crack 


2.  5  TEMPERATURE  INCREASE  OF  THE  METAL  CASE  RESULTING 
FROM  PREMATURE  EXPOSURE  OF  THE  CASE  INSULATION 

The  metal  case  beneath  the  propellant  is  protected  by  a  thin  slab  of  insulation. 

This  insulation  is  usually  only  exposed  to  the  combustion  gases  following  the 

completion  of  the  motor  firing.  However,  if  a  grain  defect  is  present  in  the 

propellant  (Figure  17),  the  case  insulation  may  be  exposed  during  the  motor 

firing,  causing  an  increase  in  the  temperature  of  the  metal  case.  The  purpose 

of  this  analysis  is  to  present  a  method  to  determine  this  temperature  as  a 

function  of  insulation  thickness  and  time  of  exposure. 

This  analysis  is  applicable  only  if  aerodynamic  heating  is  negligible.  For 
severe  aerodynamic  heating,  the  solution  of  the  governing  heat  transfer 
equations  for  the  wall  temperature  must  include  this  effect. 

However,  since  aerodynamic  heating  is  a  function  of  the  flight  trajectory, 
a  solution  to  the  problem  in  the  form  of  metal  case  temperature  as  a 
function  of  insulation  thickness  and  insulation  exposure  time  would  be 
restricted  to  a  particular  trajectory.  Therefore,  by  neglecting  aero¬ 
dynamic  heating,  a  general  first  order  estimate  of  the  effect  of  a 
particular  grain  defect  on  the  metal  wall  temperature  can  be  obtained. 

In  determining  the  temperature  of  the  metal  case  as  a  function  of  insulation 
thickness  and  insulation  exposure  time,  the  following  simplifying 
assumptions  are  employed: 

1.  Newtonian  heating  of  the  steel  case 

2.  Thermal  properties  of  the  insulation  and  steel  case  are 
invariant  with  temperature 

3.  One -dimensional  heat  transfer 

4.  Steel  case  perfectly  insulated  at  x  =  L 

5.  No  thermal  resistance  between  the  steel  case  and  insulation 

6.  Uniform  initial  temperature 

7.  Uniform  thermal  film  resistance  on  the  surface  of  the 
insulation  . 
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Figure  17.  Hypothetical  Grain  Defect 

Using  the  above  assumptions,  the  system  to  be  analyzed  is  shown  in 
Figure  18. 

The  governing  partial  differential  equation  for  heat  conduction  is 


a^T  9T 


(99) 


The  boundary  conditions  for  the  differential  equation  are: 


T  = 

T.  at  t  =  0 

A 

(100) 

T  = 

T  as  t - ►CO 

c 

(101) 

-K 

aT 

ax 

h  (T  -  T  )  at  X  =  0 
c  w 

(102) 

-K 

aT 

ax 

9T 

P  C  X  —  at  X  =  X 

^m  m  om  at  o 

(103) 

where  the  subscript  m  refers  to  the  metal  case. 
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Using  the  following  data  and  solving  for  T  at  x  =  ,  the  data  shown  on 

Figure  19  was  obtained.  The  solution,  obtained  on  the  IBM- 1620,  is  for 
material  thicknesses  corresponding  to  the  five -segment,  120-inch  motor 
with  the  following  material  properties: 


C  =  0.  35  BTU/lb  °  F 
p  =  0.  0437  Ib/cu  in. 

K  =  2.  13  X  lO'^  BTU/sec  °  F  in. 


D6AC  Steel 

C  =  0.  113  BTU/lb  °  F 
m 

p  =  0.  283  lb/ cu  in. 
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2.  6  METHOD  FOR  DETERMINING  THE  BURNING  PROPELLANT 
SURFACE  AREA  AS  A  FUNCTION  OF  TIME  BASED  ON  AN 
EXPERIMENTAL  COMBUSTION  CHAMBER  PRESSURE 

In  analyzing  a  motor  firing  in  which  the  propellant  grain  contained  numerous 

defects,  the  trace  of  the  combustion  chamber  pressure  can  be  used  to 

determine  the  burning  propellant  surface  area  as  a  function  of  time.  By 

comparing  this  burning  surface  area  against  the  burning  surface  for  the 

grain  without  defects,  the  effective  burning  surface  area  in  the  defect  can 

be  determined  approximately. 


The  system  control  volume  used  in  the  analysis  is  shown  in  Figure  20  and 
is  used  along  with  the  following  assumptions  to  obtain  the  governing 
equations  for  the  flow  model: 

A.  Assume  the  perfect  gas  equation  of  state  is  applicable 

B.  Assume  perfect  mixing  in  the  combustion  chamber 

C.  Assume  the  C  ,  C  ,  and  y  of  the  propellant  combustion 
products  in  thi  igniter  are  the  same  as  for  the  propellant 
combustion  products  in  the  main  motor  and  that  they  are 
constant 

D.  Assume  quasi-steady  flow  and  one  dimensional  isentropic 
flow  through  the  nozzle 

E.  Assume  the  heat  flux  to  the  propellant  in  the  motor  is 
negligible 

F.  Assume  the  nozzle  throat  of  the  igniter  is  choked 
initially 

G.  Assume  viscous  effects  are  negligible 

H.  Assume  the  velocity  of  the  gases  at  the  nozzle  entrance 
is  negligible  in  comparison  with  the  velocity  from  the 
nozzle  exit  . 


Due  to  the  transient  nature  of  the  ignition  process,  the  governing  equations 
which  are  used  to  define  the  flow  model  depend  on  whether  or  not  the  nozzle 
throat  of  the  main  motor  and  the  igniter  are  choked  and  whether  or  not  the 
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Figure  20.  Flow  Model 
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propellant  in  the  main  motor  is  burning.  Depending  on  the  flow  conditions, 
the  analysis  includes  the  following  regimes: 

1.  The  nozzle  throat  of  the  igniter  is  choked  but  the  nozzle 
throat  of  the  main  motor  is  not  choked,  and  the  propellant 
in  the  main  motor  is  not  burning.  The  governing  equations 
are 


w.  T. 

1  10 


w  T 
e  o 


P 

o 


(105) 


w.  =  w  +  V  p 
1  e  o 


(106) 


Aj  M 


1 


V  +  1 

2(v  -  1) 


(107) 


(108) 
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2.  The  nozzle  throat  of  both  the  igniter  and  main  motor  are 
choked,  and  the  propellant  in  the  main  motor  is  not  burning. 
The  governing  equations  are 


w.  T. 

1  10 


w  T  + 
e  o 


V 

yr 


P 

o 


w.  =  w  +  Vp 

1  e  o 


(111) 


P 

0 


RT 

o 
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3.  The  nozzle  throat  of  both  the  igniter  and  main  motor  are 
choked,  and  the  propellant  in  the  main  motor  is  burning. 
The  governing  equations  are 


V 

w.  T.  +  wT  =  wT  +  — P 
110  pc  eo  yRo 


w.  +  w  =  w+  Vp 
1  p  e  o 


(112) 


(113) 


(114) 


(115) 
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4.  The  nozzle  throat  of  the  igniter  is  not  choked  but  the  nozzle 
throat  of  the  main  motor  is  choked,  and  the  propellant  in  the 
main  motor  is  burning.  The  governing  equations  are 

V 

w,  T.  +wT  =  wT+  --;r  P 
110  pc  eo  yR  o 

w.  +  w  =  w  +  Vp 
1  p  e  '^o 


The  equations  in  each  regime  were  solved  numerically  on  the 
IBM-1620  computer. 

A  typical  solution  to  the  equations  is  shown  in  Figure  21  for  a  large 
rocket  motor.  This  motor  did  not  contain  any  grain  defects. 


Page  69 


UTC  2015-TN-l 


3.0  REFERENCES 


1.  Geckler,  R.  D. ,  "Thermal  Stresses  in  Solid  Propellant  Grain, "  Jet 
Propulsion,  February  1956,  p.  93. 

2.  Ungar,  E.  E.  ,  and  B.  W.  Shaffer,  "Thermally  Induced  Bond  Stresses 
in  Case-Bonded  Propellant  Grains, "  ARS  Journal,  April  I960,  p.  366. 

3.  Au,  N.  N.  ,  "A  Method  of  Strength  Analysis  of  Solid  Propellant  Grains," 
ARS  Paper  1062-60. 

4.  Knauss,  Wolfgang,  "Displacements  in  a  Tubular  Rocket  Grain  Under 
Axial  Acceleration,  "  GALCIT  SM  61 -10,  Guggenheim  Aero.  Lab., 

CIT,  Pasadena,  California,  April  1961 . 

5.  Pickett,  Gerald,  "Application  of  the  Fourier  Method  to  the  Solution  of 
Certain  Boundary  Problems  in  the  Theory  of  Elasticity,"  Journal  of 
Applied  Mechanics,  Vol.  11,  September  1944,  p.  A-176. 

6.  Conte,  S.  D.  ,  K.  L.  Miller,  and  C.  B.  Sensenig,  "The  Numerical 
Solution  of  Axi  symmetric  Problems  in  Elasticity,"  Proc.  of  5th 
Symp.  on  Ballis.  Missiles  and  Sp.  Tech.  ,  Vol  IV,  Academic  Press, 
1960. 

7.  Messner,  A.  M.  ,  "Propellant  Grain  Stress  Analysis, "  Bulletin  of 
17th  Meeting  JANA F-ARPA -NASA,  Vol.  II,  May  1961,  p.  149. 

8.  Timoshenko,  S.  and  J.  N.  Goodier,  Theory  of  Elasticity,  2nd  Ed.  , 
McGraw-Hill,  1951,  p.  408  . 

9.  Williams,  M.  L.  ,  P.  J.  Blatz,  R.  A.  Schapery,  "Fundamentals 
Studies  Relating  to  Systems  Analysis  of  Solid  Propellants,"  Report 
No.  GALCIT  SM61-5,  Guggenheim  Aero.  Lab.  CIT,  Pasadena, 
California,  February  1961 . 

10.  Aleck,  B.  J.  ,  "Thermal  Stresses  in  a  Rectangular  Plate  Clamped 
Along  an  Edge, "  Trans.  ASME  Vol.  71,  1949,  p.  A-118. 

11.  Courant,  R.  ,  and  D.  Hilbert,  Methods  of  Mathematical  Physics, 

Vol.  I,  Interscience  Publishers,  Inc. ,  New  York,  1955. 

12.  Lanczos,  C.  ,  The  Variational  Principles  of  Mechanics,  Vol.  4, 
Mathematical  Expositions,  Univ.  of  Toronto  Press,  Toronto, 

Ontario,  I960. 

13.  Frazer,  R.  A.,  W.  Duncan,  and  A.  R.  Collar,  Elementary  Matrices, 
Cambridge,  1952,  p.  157. 


Page  71 


UTC  2015-TN-l 


14.  Barton,  M.  V.,  Journal  of  Applied  Mechanics,  Vol.  8,  1941,  p.  A-97. 

15.  Horvay,  G.  and  J.  A.  Mirbal,  "The  End  Problem  of  Cylinders, "  Journal 
ofApplied  Mechanics,  Vol.  25,  December  1958,  p.  561  . 

16.  Kent,  C.  H.  ,  "Thermal  Stresses  in  Thin  Wall  Cylinders,"  Journal  of 
Applied  Mechanics,  APM-53-13. 

17.  Durelli,  A.  J.  and  C.  H.  Tsao,  "Determination  of  Thermal  Stresses 
in  Three-ply  Laminates,"  Journal  ofApplied  Mechanics,  Vol.  22, 

1955,  p.  190. 

18.  Zandman,  F.  ,  S.  S.  Redner,  and  D.  Post,  "Photoelastic  Coating 
Analysis  in  Thermal  Fields,"  to  be  published  ESA;  First  International 
Symposium  on  Photoelasticity,  1961  . 

19.  Williams,  M.  L.  ,  "Stress  Singularities  Resulting  from  Various  Boundary 
Conditions  in  Angular  Corners  of  Plates  in  Extension,"  Journal  ofApplied 
Mathematics,  Vol.  19,  1951,  p.  526. 

20.  Love,  A.  E.  H.  ,  A  Treatise  on  the  Mathematical  Theory  of  Elasticity, 
Dover,  1944,  p.  249. 

21.  Pickett,  S.  ,  "Fourier  Method  to  the  Solution  of  Certain  Boundary 
Problems  in  the  Theory  of  Elasticity,"  Journal  of  Applied  Mechanics, 
1944,  p.  A-176. 

21.  Kamke,  E.  ,  Differential  Gleichungen-Losungs-Methoden  und  Losungen, 
3rd  ed.  ,  Chelsea,  1959,  1.  444. 

22.  Stanford,  E.  G.  ,  and  L.  H.  Fearon,  Progress  in  Non-Destructive 
Testing,  Vol.  2,  MacMillan,  I960,  p.  39. 

23.  Welch,  J.  ,  Jr.,  "Empirical  Determination  of  Some  Mathematical 
Models  in  Viscoelasticity,"  N.  Carolina  St.  DA-01 -009-ORD-991 , 

ERD  148/1,  August  1961 .  ' 

24.  Truxal,  J.  ,  Control  Synthesis,  McGraw-Hill,  1955. 

25.  Citerley,  R.  ,  "An  Investigation  to  Simulate  the  Linear  Viscoelastic 
Behavior  of  an  Elastomer,"  MS  Thesis,  USNPS,  1961. 

26.  Grover  and  Holter,  "Solution  of  the  Transient  Heat  Conduction  Equation 
for  an  Insulated,  Infinite  Metal  Slab,"  ARS  Journal,  December  1957, 
p.  1249. 


Page  72 


UNCLASSIFIED 


UNCLASSIFIED 


